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We investigate coherent phonon propagation in a thin film of ferroelectric PbZro.2Tio.sO3 (PZT) 
by ultrafast x-ray diffraction (UXRD) experiments, which are analyzed as time-resolved reciprocal 
space mapping (RSM) in order to observe the in- and out-of-plane structural dynamics simultane- 
ously. The mosaic structure of the PZT leads to a coupling of the excited out-of-plane expansion 
to in-plane lattice dynamics on a picosecond timescale, which is not observed for out-of-plane com- 
pression. 



Oxides are attractive constituents of future nano- 
electronic devices, because of their broad spectrum of 
outstanding physical properties, such as ferroelectric- 
ity and ferromagnetism, and owing to the progress 
made in the fabrication of high quality epitaxial 
heterostructures.pQ Epitaxial strain engineering and 
the careful choice of mechanical and electrical bound- 
ary conditions enable a direct influence on these 
functionalities. [2H6] Structural defects and nanoscale in- 
homogeneities, such as dislocations and domains, typi- 
cally affect the properties of functional oxides and have 
been extensively studied by experiment and theory. [7] [8] 
Ultrafast x-ray diffraction (UXRD) emerged as a power- 
ful tool to observe lattice motion in real time[9HTT] and 
has provided a deeper insight in the structure-property 
relations of functional oxides on ultrashort timescales. 
Recent femtosecond x-ray scattering experiments on fer- 
roelectric oxides showed that electron screening induces 
an ultrafast piezoelectric response of the lattice [12] and 
that in turn the deformation leads to a change of the 
polarization. [13] However, these experiments were con- 
ducted on rather perfect epitaxial crystals. The influence 
of nano-domains has been considered in experiments on 
transient phases [14], but the role of static structural de- 
fects remained unexplored on such ultrafast timescale. 

Here we exemplify how ultrafast reciprocal space map- 
ping (URSM) using a laser-based plasma x-ray source 
yields direct additional information on the reversible in- 
plane structure dynamics in a ferroelectric perovskite 
PbZro.2Tio.sO3 (PZT) film which is solely induced by 
the existence of dislocations typical of such materials. In 
particular, the width of the PZT Bragg reflection reports 
that tensile out of plane strain leads to drastically in- 
creased damping. The energy flows into in-plane strain 
which is evidenced by the in-plane component of the re- 
ciprocal space map. Our results indicate that in mis- 
matched epitaxial films of oxide materials, with their high 
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FIG. 1. TEM and AFM images of the PZT-SRO double layer 
grown onto an STO substrate by PLD. a) The cross section 
TEM micrograph reveals that a minority of 90° a-domains 
are embedded in the matrix of the c-axis grown tetrago- 
nal PZT film, as proved by the AFM topography image in 
b) as well. Misfit dislocations (MD) and threading disloca- 
tions (TD) formed at the SRO-PZT interface and across the 
PZT film, accounting for the lateral inhomogeneity on a sub- 
100 nm length scale. 



susceptibility to the formation of domains and disloca- 
tions, in-plane phenomena emerge on a hundred picosec- 
ond timescale. URSM yields the relevant information 
on lateral lattice dynamics in such materials in which 
nanoscale inhomogeneities inherently broaden the peaks 
in reciprocal space. It is important to realize that such 
inhomogeneities are a natural paradigm in oxides orig- 
inating from competing phases with similar free energy 
rather than a result of imperfect crystal growth. [7] A bet- 
ter understanding of such time-domain effects in novel 
functional oxide materials will be important to study the 
influence of structural defects on the ultrafast response 
of collective phenomena, such as piezoelectricity. 

As a typical structure, we grew a ferroelectric layer 
of PZT and a metallic SrRuOs (SRO) electrode layer 
onto an SrTiOs (STO) substrate by pulsed laser de- 
position (PLD). [15] The transmission electron micro- 
graph (TEM) image (Fig. [l]a) shows layer thicknesses of 
dvzT = 207 nm and dsRO = 147 nm, respectively. The 
average lattice constants normal to the sample surface 
derived from static x-ray diffraction are cpzt = 4.130 A 
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and c S ro = 3.948 A. The TEM image (Fig. [T] a) fea- 
tures only few a-domains in the PZT layer, which are 
domains with a polarization vector pointing normal to 
the c-axis of the layer. [15 Accordingly, out-of-plane po- 
larized domains are called c-domains. The small amount 
of a-domains is confirmed by the very weak scattering 
observed around 3.12 A -1 in the static and transient 
rocking curves in Fig. 2]a-c). We neglect these ferroelas- 
tic domains and switching between the 90° polarizations 
states in the further discussion. Due to stress relaxation 
in the mismatched epitaxial PZT film misfit dislocations 
(MD) at the SRO interface and threading dislocations 
(TD), expanding through the complete layer, are visible 
in Fig. [T] a). [15] As a result lateral regions below 100 nm 
size are observable in PZT whereas the SRO layer is free 
of such inhomogeneities. Nevertheless, the AFM topog- 
raphy in Fig.[l]b) reveals that the mean roughness of the 
PZT surface is below 2 A. 

In order to characterize the response of the PZT film 
to ultrashort stress pulses we excited the SRO electrode 
with near infrared (800 nm) femtosecond light pulses 
(r = 40 fs) and monitored the induced lattice dynam- 
ics by UXRD experiments at a laser-driven plasma x- 
ray source (PXS) [16j [17] in a pump-probe scheme. The 
generated hard x-ray pulses (E — 8.05 keV (Cu K a ), 
r = 150 fs) were collected by a Montel multilayer mir- 
ror and focused onto the sample with a convergence of 
0.3°. The diffracted photons were accumulated with a 
CMOS hybrid-pixel area detector in classical — 20 geom- 
etry. This allowed for detecting symmetrically and asym- 
metrically diffracted x-ray photons at the same time, 
avoiding time-consuming mesh scans in order to mea- 
sure reciprocal space maps (RSM) around specific Bragg 
reflections. [TM2T] Consequently, we acquired informa- 
tion both on in- and out-of-plane structure dynamics uti- 
lizing this time-resolved version of reciprocal-space map- 
ping. The temporal overlap of the optical pump and 
x-ray probe pulses was determined in an independent 
cross-correlation experiment [22] and was set to the de- 
lay ps. 

First we discuss the conventional x-ray diffraction from 
lattice planes parallel to the surface. Fig. [2] a) shows the 
— 20 scans for three different time delays between opti- 
cal pump and x-ray probe pulses. The black line repre- 
sents the unexcited lattice and we can confirm the lattice 
constants of the three constituting materials from the 
Bragg angle. The photoinduced dynamics are evident 
from the changes of the three material specific Bragg re- 
flections, which are shown in Fig. [2]b). By fitting each 
Bragg reflection for each material with a Gaussian we 
can extract the peak width w z (t) (FWHM) and peak 
center c z (t) for each delay t in the q z dimension. For 
the later analysis it is more convenient to introduce the 
peak shift s z (t) = c z (t) — c z (0). The absorption of the 
pump pulse takes place exclusively in the SRO electrode 
layer, leading to a quasi-instantaneous temperature rise. 
The heat expansion of SRO by 0.35 % is limited by the 
speed of sound in the material and proceeds within 24 ps 
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FIG. 2. Transient x-ray diffraction measurements and static 
reciprocal space map of the PZT-SRO sample, a) — 20 
scans of the structure around the (0 2) Bragg reflections of 
PZT, SRO and STO before excitation (black) and at delays 
of maximum peak shift s z of PZT, at t = 35 ps (red) and t = 
75 ps (blue), b) Transient — 20 scans for continuous variation 
of — 10 < t < 125 ps. The horizontal dashed lines indicate 
the delays of the selected plots shown in panel a) and d). 
The solid black lines indicate the center of the Bragg peaks, 
c) The reciprocal space map of the PZT-SRO double layer 
sample before excitation features a rather broad PZT peak 
in q x direction. All peaks widths suffer from an additional 
broadening due to the convergence and energy bandwidth of 
the incident x-rays, d) Intensity of the PZT peak integrated 
over the q z dimension before excitation (black) and at delays 
of change of peak width w x of PZT, at t = 35 ps (red) and 
t = 95 ps (blue). 



evidenced by the shift to smaller Bragg angles. [28] On 
the same timescale the substrate shows a tiny shoulder 
at the high-angle side according to the compression of 
STO adjacent to the expanding metal layer. Similarly, 
the PZT film is first compressed by the strain imposed 
from the expanding SRO, however it expands after the 
strain wave is reflected from the sample surface. [23] A 
detailed evaluation and discussion will be given below. 

While this evaluation of UXRD signals is straightfor- 
ward and the example shows the power of the method, we 
now discuss how to gain the information on the in-plane 
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dynamics. The RSM before excitation is shown as a con- 
tour plot in Fig.[2]c). In general, the size of the reciprocal 
lattice points in the RSM is inversely proportional to the 
length scale of coherently scattering regions of the crys- 
tal in the according in- and out-of-plane directions. The 
additional broadening due to the instrument function of 
the x-ray diffraction setup, which is mainly given by the 
0.3° convergence and Cu K a energy bandwidth of the 
incident x-rays, can be seen in the peak profile of the 
structurally perfect STO substrate in the RSM. Fig.[2]d) 
shows the diffraction signal integrated over the q z range 
of the PZT peak. Similar to the convention above we 
define the width (FWHM) and shift for the q x dimension 
as s x (t) and w x (t). The large static value of w^ ZT is con- 
sistent with the average size of the lateral regions in the 
PZT layer of about 50 nm observed in the TEM image 
(Fig. [l]a). In crystallography this broadening of w x can 
be described by the model of mosaicitv[24]. assuming the 
crystal to consist of small mosaic blocks. These blocks 
are homogeneous in themselves but the x-rays scattered 
from different blocks do not sum up coherently. The 
in-plane size of the blocks defines the lateral correlation 
length which is inversely proportional to the broadening 
of the RSM in q x . Tilting of the blocks can give rise to 
an additional broadening. These two effects may be dis- 
tinguished by measuring a RSM around an asymmetrical 
Bragg reflection. [T8] Fig. [2] d) shows that w x increases 
considerably for snapshots recorded after the reflection 
of the strain wave at the sample surface. 

In order to discuss our experimental results we apply a 
ID model of the sample structure to simulate the lat- 
tice dynamics by a linear chain model of masses and 
springs. [25 These simulations are well established to pre- 
dict the out-of-plane dynamics but do not consider the 
in-plane dynamics directly. Therefore, we employ the 
out-of-plane phonon damping as an adjustable parame- 
ter to couple energy to in-plane motion. First, we calcu- 
late the temperature rise in the SRO after optical excita- 
tion (Fig. [3] a) from the laser fluence, absorption depth 
and heat capacity of this metal. The quasi-instantaneous 
thermal stress excites coherent acoustic phonons (strain 
waves) which are launched from the interfaces to the PZT 
layer and STO substrate. Calculated strain profiles for 
different delays are depicted in Fig. [3] b). Since heat 
diffusion from the SRO into the PZT layer can be ne- 
glected on this short time scales [26] 100 ps) the PZT 
lattice dynamics are exclusively determined by the com- 
pression wave travelling from the SRO-PZT interface to 
the PZT-air interface. Here the strain wave is reflected 
and converted into an expansion wave travelling back to 
the SRO layer and further into the STO substrate. Due 
to the good acoustic impedance matching of the three 
materials we neglect reflections at layer interfaces. 

From the simulated spatio-temporal strain map, the 
resulting transient changes of the x-ray diffraction pro- 
file in q z dimension for the PZT and SRO layers are cal- 
culated by dynamical x-ray diffraction theory. [27] Tak- 
ing the elastic constants of each material, we can use 
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FIG. 3. Simulation results for the temperature gradient and 
lattice dynamics after optical excitation, a) The optical pump 
laser causes a quasi instantaneous temperature increase only 
in the metallic SRO layer at delay zero, b) The excited lattice 
dynamics are calculated by a ID linear chain model. The 
strain profiles show a coherent soundwave with a sharp leading 
edge travelling from the PZT-SRO interface to the PZT-air 
interface. It is converted to an expansion wave with much 
stronger damping, smoothing and broadening the profile. 



the damping of coherent phonons by impurities and 
coupling to in-plane motion as adjustable parameters. 
Fig. [4]a,b,d,e) show the excellent agreement of the simu- 
lated x-ray diffraction data with the measured values for 
s z and w z . 

The peak shift s z is a measure of the average c-axis 
lattice parameter of the PZT and SRO layers, which was 
qualitatively discussed above. The change of w z essen- 
tially reflects the inhomogeneous strain, which in SRO 
is given by the short absorption length of the optical 
pump light leading to a stress exponentially decaying 
with z. Initially SRO only expands near the PZT in- 
terface. At 12 ps after excitation, the expansion wave 
has propagated through half the SRO layer, which leads 
to a maximum w^ RO (Fig. |4]d), which in fact reflects 
a splitting of the SRO Bragg peak.[25j [28] Due to the 
peak splitting the Gaussian fit indicates a compression 
of the SRO layer for sf RO , as long as only a small frac- 
tion of SRO is expanded. Similarly to SRO, w^ ZT rises 
and the peak shifts to larger angles as the compression 
travels through the PZT layer (Fig. [I]b,e). Later s^ ZT 
becomes negative when the soundwave is reflected at the 
surface (t > 50 ps) and propagates back to the substrate. 
w^ ZT is much smaller for t > 50 ps although the ten- 
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FIG. 4. a-f) Comparison of experimental data (circles) with the simulation results (solid black lines). In Panel c) and f) the 
solid lines are only guides to the eye. All grey dashed lines are the analogue to the black lines assuming constant damping 
in PZT. Panels g) and h) show the kinetic (blue), potential (red) and total (black) energy density in the according layer, j) 
Damped energy density of out-of-pane motion in PZT. 



sile strain during expansion, s^ ZT , has nearly the same 
magnitude as the preceding compression. This implies 
that the strain pulse in PZT broadens in space rendering 
the layer less inhomogeneously strained when the layer 
expands. 

In order to fit these four transient data sets for PZT 
and SRO (Fig. [4] a,b,d,e) simultaneously in our simula- 
tions we have introduced phonon damping in PZT as the 
only free parameter. The best fit for the four data sets 
was achieved when the damping factor in PZT is chosen 
50 times larger for expansive out-of-plane strain com- 
pared to compressive strain of the same magnitude. The 
according results are plotted as solid black lines in Fig. 
[|a,b,d,e). The grey dashed lines represent the results for 
a constant damping factor in PZT for expansion and com- 
pression. We can exclude pure surface scattering as the 
reason for the asymmetric damping behavior, since this 
would lead to a much smaller amplitude of s^ ZT in order 
to achieve the same decrease of w^ ZT for the expansive 
strain. The increase of the damping in PZT is visualized 
in Fig. [3] b) where the coloured lines (t > 50 ps) show a 
smooth out-of-plane expansion, whereas the greyish lines 
(t < 50 ps) show the inhomogeneous compression before 
the reflection at the surface changes the sign of the strain 
wave. 

From our lattice dynamics simulations we can also 
determine the kinetic (blue), potential (red) and to- 
tal (black) energy density of the out-of-plane coherent 



phonons in each layer (Fig. [4] g-h). We introduced 
the phonon damping in PZT to couple energy to lat- 
eral phonons. This energy is essentially the difference 
of the total energy of out-of-plane coherent phonons in 
PZT with and without damping. The result is plotted in 
Fig. |4]j) where the grey dashed line corresponds again 
to the case of constant damping. The increase of the lat- 
eral energv density in PZT (Fig. [|j) goes along with a 
considerable increase of w^ ZT which reflects a change of 
the inhomogeneity in plane, probably because the lateral 
blocks develop an inhomogeneous in-plane strain that is 
dynamically coupled to the out-of-plane motion accord- 
ing to the Poisson ratio (Fig. |4]f). We do not observe 
the converse effect during the compression of the PZT 
layer (t < 50 ps). The horizontal arrow in Fig. [4]f) in- 
dicates the timescale of the build-up of the lateral strain 
of approx. 50 ps. We can link this timescale to a lateral 
lengthscale of approx. 200 nm via the sound velocity of 
PZT of 4.6 nm/ps. This timescale agrees well with the 
in-plane block size observed in the TEM image in Fig. 
[l]a) and w^ ZT of the static RSM in Fig. [|c). 

The analysis of the measured signal alone already sug- 
gests the following interpretation: The expansion of SRO 
sends a compression wave into PZT. The in-plane mosaic- 
ity/nanoinhomogeneity is unchanged during this period. 
When the strain changes sign upon reflection at the sur- 
face PZT expands and according to Poisson ratio, the 
mosaic blocks must now laterally contract. The in-plane 
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inhomogeneity is increased as millions of in-plane con- 
traction waves start at all the lateral dislocations. We 
conclude that only out-of-plane expansion of PZT cou- 
ples energy to in-plane dynamics, and that this effect is 
essentially suppressed for out-of-plane compression, since 
this would have to expand the mosaic blocks, which is 
sterically forbidden by the adjacent blocks. This com- 
pares favourably with our simulations of the out-of-plane 
lattice dynamics, requiring an increased damping for the 
expansion wave in the PZT, which can be understood as 
an increase of the coupling between in- and out-of plane 
lattice motions. 

In conclusion we have demonstrated the first measure- 
ment of the lattice dynamics in a structurally imperfect 
thin film by ultrafast reciprocal space mapping (URSM). 



We not only extract the changes of the lattice constants, 
i.e. the expansion and compression of materials perpen- 
dicular to the surface. In addition, we quantify the cou- 
pled response in plane, which turns out to be significantly 
enhanced for out-of-plan expansion, as it provides the in- 
plane contraction necessary for the atoms to start mov- 
ing. URSM will be an important method for understand- 
ing the ultrafast response of oxide crystals with their nat- 
ural tendency to form nanoscale inhomogeneities. 
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